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Temperature and abundance profiles of hot gas in galaxy groups - 
I. Results and statistical analysis 
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ABSTRACT 

The distribution of metals in groups of galaxies holds important information about the chemi- 
cal enrichment history of the Universe. Here we present radial profiles of temperature and the 
abundance of iron and silicon of the hot intragroup medium for a sample of 15 nearby groups 
of galaxies observed by Chandra, selected for their regular X-ray morphology. All but one 
group display a cool core, the size of which is found to correlate with the mean temperature 
of the group derived outside this core. When scaled to this mean temperature, the temperature 
profiles are remarkably similar, being analogous to those of more massive clusters at large 
radii but significantly flatter inwards of the temperature peak. The Fe abundance generally 
shows a central excess followed by a radial decline, reaching a typical value of 0.1 Zq within 
rsoo, a factor of two lower than corresponding results for clusters. Si shows less systematic 
radial variation, on average displaying a less pronounced decline than Fe and showing evi- 
dence for a flattening at large radii. Off-centre abundance peaks are seen both for Fe and Si 
in a number of groups with well-resolved cores. Derived abundance ratios indicate that su- 
pernovae type la are responsible for 80 per cent of the Fe in the group core, but the type II 
contribution increases with radius and completely dominates at 7-500 • We present fitting for- 
mulae for the radial dependence of temperature and abundances, to facilitate comparison to 
results of numerical simulations of group formation and evolution. In a companion paper, we 
discuss the implications of these results for feedback and enrichment in galaxy groups. 
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1 INTRODUCTION 

Hot X-ray emitting gas constitutes the dominant baryonic com- 
ponent in groups and clusters of galaxies. The metals in this 
intracluster medium (ICM) are believed to originate mainly 
in mat erial ejected from g roup and cluster galaxies by super- 
novae jAmaud et al.l 1 19921) . with a smaller portion driven out 
by galaxies through ga laxy-galaxy a nd ggJaxy-ICM in teractions 
jPomainko et alj I2OO6I, but see also IMoII et alj l2007h . Besides 
galactic stellar winds, other non-gravitational processes have also 
played a role in redistributing the enriched gas, and more generally, 
in shaping the thermodynamic properties of the ICM. These pro- 
cesse s include radiative cooling ( Sanders on. Poimian & O'SullivanI 
l2006h . and he ating and mixing of ga s by a c tive galactic nu- 



clei Jcros ton. Hardcastl e & Birkinshaw 20051: McNamara et al.l 



20051: lO' S ullivan. Vrtile k & Kempnej|2005l : iRebusco et al.ll2006l : 



Roediger et alj |2007) . The spatial distribution of metals, and to 



some extent the temperature structure of the ICM, reflect the action 
of such processes, thus offering insight into the importance and na- 
ture of the mechanisms, other than gravity, that have established the 
properties of the ICM at the present epoch. Groups of galaxies are 



E-mail: jesper@star.sr.bham.ac.uk 



particularly interesting in this context because such processes are 
expected to be relatively more important in low-mass systems due 
to their shallower gravitational potentials. 

Iron and silicon comprise the major diagnostic elements for 
studies of the ICM abundance distribution, as these elements give 
rise to the most prominent lines in the soft X-ray spectrum of ther- 
mal plasmas, and are among the elements for which there is a 
distinctively different yield in type la and type II supernova ex- 
plosions according to st andard supernova (SN) models (see, e.g., 
'Baumgartner et al.'2005 and references therein). Whilst iron is pre- 
dominantly produced by SN la, silicon is more evenly mixed be- 
tween the two SN types. The ratio of Si to Fe abundance therefore 
provides valuable information on the relative importance of SN la 
vs. SN II in enriching the ICM. 

Abundance and temperature pro files of massive clusters have 
received c onsidera b le atte n tion (e.g..lDe Grandi & Molendill200ll: 
Baumgart ner et alj l2005l: Ivikhlinin et alj l2005l : IPiffaretti et al.1 
120051 : iDonahue et alj|2006l) . It is now well established that clus- 
ters with cool ce ntral regions typically exhibit high central Fe 
abundances (e.g., De Grandi & Molendil I2OOII : llrwin & BregmanI 
I2OOIL iDe Grandi et al J |2004) which can be attributed predomi- 
nantly to SN la enrichment by the stellar population of the cen- 
tral bright cluster galaxy jFinoguenov. David & PonmanI bOOd : 
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iBohringer et"ai]|2004 lOe Grandi et"ai]|2004h . Clusters with no ev- 
idence for strong central cooling tend to display a more uniform Fe 
distribution. 

At large radii (r > 0.5r2oo, where r2oo is the radius en- 
closing a mean density of 200 times t he critical density), Zpe 
typically drops to a value ~ 0.2 i Finoguenov et al.' '200dl; 



De Grandi & Molendil I2OO1I : iTamura et a l. 2004; De Grandi et al. 
2004) . As suggested by. for example, Finogue nov et al.H2000f) . this 



could result from SN Il-dominated enrichment at an early stage in 
the cluster formation via energetic starburst winds. Such a picture 
receives support from the observation that a sub stantial fraction of 
the ICM metals wer e already in place by 2 ~ 1 jTozzi et al.ll2003l : 
iBalestra et al]|2007h . 

The situation in lower-mass systems is much less clear, de- 
spite their relatively larger importance for the cosmic baryon 
and metal budgets. Historically, this has reflected the diffi- 
culty in assembling large samples of groups bright enough 
to allow detailed, spatially resolved X-ray spectroscopy. Ear- 
lier works based on ROSAT and ASCA data typically featured 
only a handful of kT < 2 keV groups, oft en included as 
part of larger cluster-dominat e d sam p les (e.g 
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19981: iFinoguenov & Ponman 1 19991 : iFinoguenov e t al. 2000; 
Finoguenov. Arnaud & David 2001al) . Such st udies suggested 



a Zs i /Z-pe ratio roughly solar in groups jFukazawa et al 
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1999I). about half the value found 

jMushotzkvetalJI 19961) . with an indication that Si/Fe rises above 
the solar ratio outside group cores (e.g. lFukazawaetai]|l998l : 
IFinoguenov & PormianI [l993) . However, the sensitivity and spa- 
tial resolution of ROSAT and ASCA imposed natural limits to the 
amount of detail one could obtain on the spatial variation of Fe and 
Si and hence the relati ve importance of the tw o supernova types 
at different group radii ( IFinoguenov et al.ll2002h . For the same rea- 
sons, it is possible that the so-called Fe and Si biases i lBuotel2000i]) 
may have affected some of these earlier results, at least for groups 
with thermally complex core regions. 

With the current generation of X-ray telescopes, the situa- 
tion has improved dramatically. Recent XMM observations have 
shown that temperature profiles of cool-core groups display a 
striking similarity and typically peak at a radius r ~ O.OSrioo 
jGastaldello et alj|2005l) . which corresponds to r ~ O.lrsoo for 
the mass profiles derived for these groups. An other recent study, 
also b ased on XMM observations, is the work of IFinoguenov et al] 
( l2007l) . who investigated the distribution of temperature, density, 
metal abundance, and derived quantities using 2-D spectral map- 
ping of a sample of 14 groups. This work showed that, as for clus- 
ters, central metal excesses are present in groups with a cool core, 
a result already hinted at in earlier ASCA studies teuotdl20 00al. 

While Chandra lacks the sensitivity of XMM-Newton, it has 
a superior spatial resolution, enabling a more robust removal of 
point sources, and also experiences a lower and more stable in- 
strumental background. Chandra data can therefore complement 
XMM observations of groups, particularly in the group cores, but 
can also serve as an important validity check on XMM results for 
the X-ray faint outskirts of groups due to the lower background. 
Here we utilize this to derive radial profiles of temperature and the 
abundance of Fe and Si from Chandra observations of a sample 
of 15 groups. Particular emphasis is given to combining the results 
obtained for individual groups to explore statistical trends within 
the sample. One purpose of doing so is to obtain a clearer picture 
of the content and behaviour of metals in group outskirts, where 
most of the intragroup gas resides but where the low X-ray surface 



brightness often prohibits robust constraints for individual systems. 
In addition, this approach allows an investigation into the signa- 
tures of feedback in group gas, including the origin of metals in 
SN and the redistribution of gas by active galactic nuclei (AGN), 
in more detail than is possible on a group-by-group basis. Further- 
more, many numerical simulations of groups and clusters now in- 
clude chemical evolutio n of the ICM through t he actions of galac- 
tic starburst winds (e.g., Tornatore et al. 2004; Romeo et al. 200^; 
ICora 2006 ), AGN outflows (Moll et al . 2007), and galaxy-ICM in- 
teractions toomainko et al.ll2006l) . Therefore, it also seems timely 
to provide a firm observational baseline against which to test the 
input physics in such simulations, from the perspective of their pre- 
dictions for the chemical (and thermal) properties of the ICM in the 
group regime. 

The paper is structured as follows. In Section|2|we describe the 
group sample and data analysis procedure, and SectionjSjpresents 
the main results for each group. Section |4| discusses the impact of 
systematic errors on our results. In Section|5]we stack the results 
for individual groups, to investigate and quantify systematic trends 
across the sample. Section |6| summarizes the results and presents 
our conclusions. In a companion paper (Rasmussen & Ponman, in 
prep.; hereafter Paper II), we discuss the implications of these re- 
sults for the history and nature of galactic feedback and chemical 
enrichment within the sample. Ho = 70 km s~^ Mpc~^ is assumed 
throughout, and all errors are given at the 68 per cent confidence 
level unless otherwise stated. 



2 GROUP SAMPLE AND DATA ANALYSIS 

2.1 Sample selection and data preparation 

Our sample is predo minantly based on the 25 GEMS groups 
from the 'G' -sample of Osmond & P orrniarj j2004h . This subsam- 
ple comprises those GEMS groups that display group-scale X-ray 
emission (spatial extent of r > 60 kpc), distinguished this way 
from systems that remain X-ray undetected and those in which the 
diffuse emission can be associated with the central galaxy rather 
than with the group as a whole. From these 25 groups, we effec- 
tively imposed a flux limit on the sample by only selecting systems 
with available Chandra archival data containing more than 6,000 
ACIS net counts from diffuse emission. This was done in order to 
enable robust, spatially resolved spectroscopy in a number of re- 
gions for each group. In addition, only groups more distant than 
20 Mpc were included, to ensure that data extended well outside 
the central group galaxy. As one of our aims is to investigate signa- 
tures of galactic feedback in group gas, and since we wish to study 
the radial distributions of gas properties, this study requires systems 
which appear reasonably undisturbed by recent merger activity and 
associated gas mixing. We therefore also excluded NGC 5 17 1 from 
th e sample, an obviously unrelaxed system already studied in detail 
bv lOsmond, Ponman & Finoguenov! j2004h . 

This left II groups, to which we have added a further four 
groups with data in the Chandra archive that conform to the 
same criteria. The resulting sample includes two Hickson compact 
groups (HCG 42 and 62) and a foss il group (NGC 741) according to 
the definition of Ijones et and spans a temperature range 

kT — 0.3-2.1 keV, with a median temperature of T = 1.06 keV 
and median redshift of z — 0.014. It is worth emphasizing that the 
selected groups are all reasonably X-ray bright and that the sam- 
ple is by no means statistically representative. In fact, the invoked 
net counts criterion alone eliminated 13 of the 25 GEMS systems 
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from which our sample was largely drawn, resulting in a sample 
lower limit of Lx ~ 5 x 10*^ erg s~^. Three of our 15 groups 
are also included in the XMM sample of Finoguenov et al. ( 2006), 
with a further 6 included in the XMM sample of Finoguenov et al. 
( l2007h . enabling cross-comparison of the results obtained from the 
two instruments. 

For all our data sets, calibrated event lists were regenerated 
using CIAO v3.3. For Very Faint mode observations, the standard 
additional background screening was carried out. Bad pixels were 
screened out using the bad pixel map provided by the pipeline, 
and remaining events were grade filtered, excluding ASCA grades 
1, 5, and 7. Periods of high background were filtered using 3- 
a clipping of full-chip lightcurves, binned in time bins of length 
259. 8-s and extracted in off-source regions in the 2.5-7 keV band 
for back-illuminated chips and 0.3-12 keV for front-illuminated 
chips. Blank-sky background data from the calibration database 
were screened and filtered as for source data, and reprojected to 
match the aspect solution of the latter. 

Point source searches were carried out with the CIAO task 
'WAVDETECT' using a range of scales and detection thresholds, 
and results were combined. Source extents were quantified using 
the 4(7 detection ellipses from 'WAVDETECT', and these regions 
were masked out in all subsequent analysis. Diffuse emission from 
individual group galaxies other than the central early-type present 
in all groups was also masked out. 

To provide a rough idea about the spatial distribution of 
X-ray emission in each system, images of all groups are pre- 
sented in Fig. [T] showing the adaptively smoothed, background- 
subtracted and exposure-corrected 0.7-2.5 keV emission on the 
relevant central CCD (the ACIS-I array for NGC 6338 and 
NGC 7619, ACIS-S2 for NGC 2300, and S3 for all others). 
T hese images were generated following t he procedure outlined 
in lRasmussen. Ponman & Mulchae\j j2006l) . Although small-scale 
substructure is visible in some of the systems, with the second 
brightest group galaxy also seen in X-rays in a few cases, most 
of the groups are evidently reasonably dynamically relaxed, sup- 
porting the use of one-dimensional radial profiles to describe the 
spatial variation in hot gas properties of the systems. 

2.2 Spectral fitting procedure 

All regions used for spectral extraction were centred on the peak of 
the diffuse X-ray emission, which in most cases coincided to within 
1-2 arcsec (and always within 5 arcsec) with the optical position of 
the brightest group galaxy as listed in the Hyperleda database. For 
each CCD in each data set, the appropriate blank-sky data were 
scaled to match off-source count rates in the 10-12 keV band and 
then used to generate background spectra. All spectra were fitted 
in XSPEC vll.3 in the 0.7-5 keV band, with the low-energy cut- 
off chosen such as to reduce contamination from residual Galactic 
soft X-ray emission not accounted for by our blank-sky background 
subtraction (see Section |4] for details). Indeed, no systematic soft 
residuals between source and background spectra were seen within 
this energy band for any of the groups. 

To put all groups on an equal footing, all radii were normalized 
to rsoo, using 

r50o(r) = 0A5T^^^ h^^ Mpc, (1) 

which is an empirical relation based on data outside the cluster core 
for a sample of 39 groups and clusters, including eight syste ms with 
r < 1.5 keV jFinoguenov, Reiprich & B6hringeijr2001bh . In or- 
der to obtain a representative mean temperature (T) for use with 



equation (TJ, as well as a corresponding mean abundance l^Z), and 
to create spectral weights for subsequent fits in individual annuli, 
a source spectrum was first extracted within the radial range 0. 1- 
0.3 rsoo. This range was chosen because it extends to the maximum 
radius covered by ACIS for the more nearby systems, while in all 
cases being well outside the optical extent of the central galaxy. 
The latter was quantified as -D25, the ellipse outlining a _B-band 
isophotal level of 25 mag arcsec"^ . The associated fit results should 
therefore remain largely uncontaminated by X-ray emission from 
the interstellar medium and any unresolved point sources in the 
central galaxy. As we will see, this choice of radial range also ex- 
cludes any cool core, in accord with our usage of equation (TJ. The 
fits to these spectra all adopted an absorbed thermal APEC model 
in XSP EC, assuming the abundance table from .Grevesse & Sauvall 
( Il998l) . Using the derived value of (T) with equation ([T), rsoo was 
then re-evaluated iteratively, until the resulting values of (T) var- 
ied by less than 1 per cent. Throughout, we will refer to the values 
of (r) and {Z) obtained in this way as the mean temperature and 
abundance of each group. 

Spectra were then extracted in concentric annuli containing 
at least 2000 net counts in the 0.7-1.5 keV band, and accumu- 
lated in bins of at least 25 net counts. Spectral response products 
were weighted by a model of the 'mean' spectrum extracted above. 
Again, an absorbed APEC model was first fitted to the data, yield- 
ing a mean temperature and abundance for the relevant annulus. A 
second fit was then performed using an absorbed VAPEC model, 
with all element abundances tied together except for Fe and Si. For 
spectra extracted inside the region covered by the central early-type 
galaxy, a power-law of fixed index F = 1.72 was added to the 
spectral model, to account for any unresolved contribution from 
low-mass X-ray binaries. 

To investigate the presence of multiple thermal components 
and add ress the impor tance of systematic errors associated with the 
Fe bias ( lBuotel2000ah . we also fitted a 2-T plasma model to all bins 
inside the optical extent of the central galaxy, with both the Fe and 
Si abundance tied for the two thermal components. Where either 
the resulting Fe or Si abundance was significantly different from the 
\-T result (at 90 per cent confidence), we adopted the best-fitting 
2-T values for both metals, confirming that a 2-T model provided a 
better fit in all cases. As discussed in more detail in Section|4l this 
was relevant for only five groups. We note that the adoption of the 
2-T results for the abundance in these five cases does not affect (Z) 
significantly, as this quantity has been derived within 0.1-0.3 rsoo, 
which is generally well outside the region for which 2-T models 
provide superior fits in the relevant groups. 

Owing to the layout of the ACIS chips, we do not have full 
azimuthal coverage of the groups at large angular distances from 
the group centres. Our results at large radii could therefore be sen- 
sitive to azimuthal variations in fitted parameters. In an attempt to 
accommodate the systematic uncertainties associated with this, we 
used azimuthal variations at smaller radii to gauge the variations 
further out. For each group, a number of spectra were extracted 
from an annulus outside D25 extending to the maximum radius of 
full coverage. To mimic the situation at large radii, only data within 
a narrow strip across this annulus were used, masking out the re- 
maining regions to achieve, if possible, 2000 net counts within the 
strip, as for the full-coverage annuli. This was repeated a number of 
times with varying strip position angles. The resulting spectra were 
then each fitted with a VAPEC model as described above, and the 
resulting standard variations on relevant parameters (T, Z^c, Zsi) 
adopted as the systematic uncertainties at large radii. The mean un- 
certainties resulting from this approach were 8, 34, and 58 per cent 
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Figure 1. Adaptively smoothed 0.7-2.5 keV images of the central region of all groups. Images of NGC 6338 and NGC 7619 are 14.4 arcmin across, all others 
7.2 arcmin. Dashed ellipses outline the D25 ellipse of the brightest group galaxy. 



on jT, Zfc, and Zsi, respectively. These errors have been added in 
quadrature to the statistical errors obtained for all annuli with less 
than 75 per cent azimuthal coverage (for most groups just one or 
two radial bins). 

Table[T]records the position of the X-ray peak of each group, 
group distance D corrected for Virgocentric infall, the useful ACIS 
exposure time tcxp, {T) and (Z) from the 'mean' spectrum, rsoo 
from equation {T](, and the Galactic value of the absorbing column 
density A^h from lOickev & LockmanI jl99Ct) as adopted in the fits. 
Note that the statistical uncertainty on rsoo based on that of (T) 
alone is typically on the order of just a few per cent. The actual 
uncertainties on rsoo are therefore likely to be dominated by sys- 



tematics associated with the use of equation ([TJ, but a quantification 
of this is beyond the scope of this paper. 



3 RESULTS 

In this section we present the results for the projected radial profiles 
of temperature, Fe and Si abundance, and the ratio Zsi/Zpc- As one 
of our aims is a statistical investigation of radial trends across the 
sample, the profiles have not been deprojected, since this would 
remove the statistical independence of the individual radial bins. 
However, we discuss the impact of deprojecting the data in Sec- 
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Table 1. Properties of the group sample. 



Group 



RA 
(J2000) 



Dec 
(J2000) 



D 
(Mpc) 



t'CXp 

(ks) 



(T> 
(keV) 



{Z) 
(Z©) 



»"500 
(kpc) 



(10^0 cn 



NGC 383 


01 07 24.9 


+32 24 45.8 


73 


42.5 


1 


NGC 507 


01 23 39.2 


+33 15 19.6 


69 


14.3 


1 


NGC 533 


01 25 31.4 


+01 45 33.7 


76 


33.6 


1 


NGC 741 


01 56 21.0 


+05 37 42.6 


79 


28.4 


1 


NGC 1407 


03 40 11.8 


"18 34 48.9 


26 


37.7 


1 


NGC 2300 


07 32 19.2 


+85 42 31.3 


30 


43.9 





HCG42 


10 00 14.2 


-19 38 12.5 


64 


30.5 





MKW4 


12 04 26.9 


+01 53 44.9 


96 


29.3 


1 


NGC 4125 


12 08 05.7 


+65 10 27.7 


22 


60.0 





NGC 4325 


12 23 06.7 


+10 37 16.0 


117 


27.9 





HCG 62 


12 53 05.8 


-09 12 16.0 


74 


47.1 


1 


NGC 5044 


13 15 24.0 


-16 23 06.4 


33 


18.0 


1 


NGC 5846 


15 06 29.4 


+01 36 23.3 


30 


22.8 





NGC 6338 


17 15 22.9 


+57 24 40.5 


127 


33.5 


2 


NGC 7619 


23 20 14.5 


+08 12 24.9 


52 


25.6 


1 



r + 0.04 

^-0.06 

r. + 0.03 

^-0.03 
,+0.05 
^-0.05 
,+0.14 
^-0.12 
,+0.07 
^-0.09 
B+0.04 
'-0.03 
„+0.05 
^-0.05 
a+0.07 
'-0.09 
O+0.12 
^-0.05 
□+0.02 
='-0.02 

r, + 0.03 

■^-0.03 
,+0.03 
^-0.03 
K+0.04 
"-0.03 
,+0.19 
^-0.07 
K+0.07 
"-0.03 



39+0 0'' 
'^■•^='-0.07 
-1+0.07 
-0.06 



0.40^ 



28+0-07 



0.16 
0.15 
0.17; 
0.25 
0.49 
0.19 
0.39 



+0.08 
0.06 
+0.09 
0.06 
+0.07 
0.05 
+0.31 
0.10 
+0.07 
0.07 
+0.17 
0.07 
+0.08 



-0.06 
:,+0.02 
-0.03 
-+0.04 
-0.03 
16+006 

-+0.08 
-0.05 

r, ,0 + 0.05 

^■^'-'-0.05 



0.12: 
0.25^ 



0.25^ 



578 


5.41 


513 


5.23 


497 


3.07 


536 


4.44 


452 


5.42 


397 


5.52 


402 


4.78 


600 


1.88 


259 


1.82 


448 


2.14 


450 


3.06 


476 


4.94 


366 


4.24 


657 


2.60 


463 


5.04 



tion|4] showing that this would not compromise our general con- 
clusions. 

In order to compare observed abundances and their ratios to 
those expected from SN la and SN II enrichment, we adopted 
SN model yields taken from the literature. For SN la, the cur- 
rently favoured dela yed-detonation WDD2 scenario is assumed 
jlwamoto et alj[l999l) . while for core-collapse supe movae (includ- 
ing typ es lb and Ic), we have used the recent yields of lNomoto et alj 
( l2006h . The latter were mass-averaged using a Salpeter initial mass 
function (IMF) for the progenitors over the range 10-50 M© ( see 
iFinoguenov, Burkert & Bohringedl2003l and lde Plaa et al ] |2007l for 
a justification of this choice of IMF), so that 



Mi 



(2) 



is the mass of the i-th element produced in a star of mass m, 
with X = 1.35. We have assumed a fractional metal mass of 
the SN II progenitors of Z = 0.004 (i.e. Z « 0.2 Z©), but 
the exact choice has no strong beari ng on the resulting Fe and Si 
yields (see, e.g.. Ide Plaa et al.l l 20071) . The resulting average stel- 
lar yields are Mpe « 0.79 M© and A'/si « 0.21 Mq for Fe 
and Si from SN la, and Mpo ~ 0.08 Mq and Msi « 0.12 Mq 
for core-collapse supernovae including SN II. We note that the 
yield uncertainties for SN II are substantial, a factor ~ 2 for Fe 
and ~ 30 per cent for Si; the yields adopted here are bracketed 
by the results of the different mode l s cons idered by, for example, 
iGibson. Loewenstein & Mushotzkvl ( 1 19971) . and so can be viewed 
as representative of typical values in the literature. With the adopted 
abundance table, the nominal abundance ratios expected for pure 
SN la and SN II enrichment are Zsi/Z^c ~ 0.46 Zsi.o/Zpo,© 
(SN la) and 2.62 Zsi.o/Zpc,© (SN II). Note that this implies SN II 
dominance (in absolute numbers) for Zsi/^Fc -Sj 0.65, and that a 
solar abundance ratio requires 3.4 SN II per SN la. 

By adopting the SN la yields o f llwamoto et"ai] ( ll999l) . we re- 
strict ourselves to the 'standard' picture of a reasonably homoge- 
neous population of SN la for simplicity. We note, however, that 
recent developments indicate that the situation could be more com- 
plex, with evidence for a diversity in SN la yields in cluster cores 



('Finoguenov et al."20 02|) and for two d istinct populations of SN la 
(iMannucci. Delia Valle & Panagiall2006) . 



3.1 Temperature profiles 

Radial temperature profiles for all groups derived from single- 
temperature models are presented in Fig. (2] with the radii of all 
annuli normalized to rsoo. It is worth noting that emission is de- 
tected to fairly large radii in most of the groups, with the midpoint 
of our outermost radial bin corresponding to roughly 75 per cent of 
rsoo . This can be ascribed to the relatively low and stable Chandra 
background. 

In general, the temperature profiles appear remarkably similar, 
with all but one system (NGC 4125) showing evidence for cooler 
gas in the central regions. The temperature profiles typically peak 
around r ~ O.lrsoo and decline at large radii. There is also an 
indication that the profiles flatten in the centr al regions, in gen - 
eral agremeent with the group results of Finog uenov et al.l yOOj). 
Compared to analogous profiles for relaxed cluster s, this contrasts 
with the projected Chandra temperat ure profile s oflVikhlinin et al] 
12005), the deprojected ones of Sande rson et al.li 2006h . and the de- 
projected XMM profiles of Piffaretti et al. "(2005)7 for which there 
is no clear evidence for a flattening at small radii. Defining the gas 
en tropy as S = TnZ^^ ^ , the deprojected cluster entropy profiles 
of iDonahue et al. I (I2OO6I) do show a central flattening, but this was 
attributed to a central core in the density distribution rather than a 
flattening of the temperature profile. 

The typical location of the temperature peak s in our sample is 
in exc ellent agreement with the XMM results of iGastaldello et al] 
( |2005|) . who also found temperature peaks at r ~ O.lrsoo for 
their cool-core groups and poor clusters. Their sample has four 
groups in common with ours (NGC 533, MKW 4, NGC 4325, 
and NGC 5044), and has rsoo evaluated directly from the mass 
profiles derived for each group. The good agreement with our re- 
sults therefore also suggests that our estimates of rsoo are gener- 
ally not seriously biased. As also noted by Gastal dello et alj j2005h 
for their sample, the temperatures peak at larger overdensity than 
in clusters, where the profiles turn over at r ~ 0.15r5oo ~ 



he pr 

O.lr2oo jPiffaretti et alj]2005l ; IVikhlinin et al] |2005| ; IZhang et al.1 
I2OO6I : ISanderson et alJbom It is clear, however, even when nor- 
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Figure 2. Projected radial temperature profiles of all groups. Vertical dotted lines mark the seini-major axis of the D25 ellipse of the central group galaxy. 



malized to rsoo, that the exact radius r(rmax) at which the peak 
occurs (in other words, the radial extent of the cool core) varies 
slightly from group to group. Naturally, this result will to some ex- 
tent depend on the spatial resolution resulting from our binning cri- 
terion and thus the surface brightness of the group emission, but we 
note that the adopted spatial resolution does not depend systemat- 



ically on, for example, group temperature. We discuss the location 
of the temperature peak in more detail in Section lSTI 
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3.2 Abundance profiles 

The abundance profiles of Fe and Si are plotted in Fig. [3] In all 
systems except NGC 2300, and possibly the non-cool core group 
NGC 4125, there is clear evidence for a central excess of iron, 
as also seen for cool-core clusters. This central excess often ex- 
tends well beyond the optical extent of the central galaxy, and 
in some cases shows an off-centre peak inside -D25, with a de- 
cline in abundance in the very centre. NGC 533 is perhaps the 
clearest example of this, but the majority of groups with a well- 
resolved core (roughly 8/11) shows some evidence of this fea- 
ture. This Fe peak is typically accompanied by a similar fea- 
ture in the Si distribution, most clearly in HCG 42, MKW 4, 
and HCG 62, with a hint of similar behaviour in a few more 
groups. Potentially, the peak could be ascribed to Fe/Si bias in- 
side the central galaxy, lowering the derived abundance in the very 
core, but as discussed below, we find little evidence that this ef- 
fect is important for our analysis in general. In addition, such a 
peak is commonly seen also in well-resolved cluster cores where 
the Fe bias is generally expect ed to be l ess important, e.g., Cen- 
tauru s jSanders & Fabianl2002h. Pe rseus jSanders. Fabian & Dunn 
2005t) . AI795 jEttori et all I2OO2I) , and A2I99 j Johnstone etal 
2002 h. Alternative explanatio ns for the p eak include A GN redis- 
tribut ion of enriched gas (e.g.. lMathews et al. 2003; Reb usco et al.l 
|2006|). resonant scattering of emission line photons in group cores 
jGilfanov. Sunvaev & Churazovlll987[). or gravitational s edimenta- 
tion of helium nuclei in the core tettori & FabiM]|2006l) . We will 
discuss these possibilities in detail in Paper II. 

At large radii, a common feature of the Fe distributions is a 
drop in Z^c to a value of ~ 0.1 Z©. In some groups, a rise in Zsi 
is also seen, largely independent of the behaviour of the Fe profile. 
We note that there is no indication that derived profiles for the two 
Hickson groups or the fossil group NGC 741 differ significantly 
from those of the rest of the sample in this context. 

The silicon-to-iron ratio Zsi/Z-pc is plotted for all groups in 
Fig. |4] In many cases the ratio in the central bin is seen to be con- 
sistent with the Solar ratio, with a contribution from SN II formally 
required for about half the groups. For the other half, there is clear 
evidence for subsolar abundance ratios, consistent with enrichment 
by SN la alone. We do not believe there is any fundamental dif- 
ference between the two subsamples if split according to this, as 
Zsi/Z-pc in the central bin clearly spans a continuous range of val- 
ues rather than showing a bimodal distribution. Given the varia- 
tion in spatial resolution among the groups and the uncertainties 
on Zsi/Z-pc, we caution further against such an interpretation at 
this stage. For example, we find no evidence that the two subsam- 
ples cluster differently in the parameter space spanned by (T), {Z), 
and if-band luminosity of the central galaxy. Nevertheless, we will 
revisit this issue in Paper II on the basis of integrated Fe and Si 
masses inside _D25 . 

As discussed by iFinoguenov et alj ^20021) . there is observa- 
tional evidence for diversity in SN la yields in the cool cores of 
more massive clusters, with solar Si/Fe ratios arising entirely from 
SN la. While our results allow for this possibility, we cannot eas- 
ily confirm that this result extends down into the group regime, 
as it is based on the abundances of additional elements such as 
S, Ca, and Ar. At the lower temperatures prevalent in groups, the 
abundances of these elements, with prominent emission lines in the 
range E ~ 2.5 — 5 keV, are not in general robustly constrained. 
In the outer regions, the interpretation of our data remains unaf- 
fected by this; here a contribution from SN II is universally required 
within our sample, and the abundance ratios, given the typical un- 



certainties in SN II yields of a factor ~ 2 (e.g. lGibson et al.lll997h . 
are consistent with pure SN II enrichment in all cases. 

3.3 Notes on individual systems 

Many of the groups exhibit properties which warrant some 
additional comments. Below, we briefly discuss each group in the 
sample in turn, with emphasis on the derived abundances. Where 
relevant, we have also attempted a comparison of the derived 
Fe profile to existing measu r ement s , in p articular the 2-D XMM 
results of 'Finogueno v et al.l feOOd l2007h (F06 and F07 in the 
following), while correcting for differences in adopted abundance 
tables. It should be noted, however, that due to our incomplete 
azimuthal coverage at large radii, direct comparison is not always 
straightforward for sources with some degree of asymmetry in the 
Fe distribution. 

NGC 383: The central galaxy hosts the bright radio source 
3C 31, representing one of the few examples of a central group 
galaxy showing a clearly id entified large- scale radio jet with a 
well-constrained jet power dLaing & Bridle! 1200 2). Despite this, 
the T- and Z-profiles of this group do not exhibit any clear 
distinguishing features compared to the rest of our sample. For the 
temperature profile, this result is in accord with the group results 
of Jetha et al. (2007), which show that the slope of T(r) inside the 
cool core does not show any strong dependence on the radio power 
of the central radio source. 

NGC 507: An example of a group in which radio lobes are present 
around the central galaxy (Kraft etal. 2004) and in which the 
highest-metallicity gas is displaced from the group core, perhaps 
as a consequence of radio outbursts lifting enriched gas out of the 
core. If so, one might expect to see some degree of asymmetry in 
the Fe distribution in the core. The 2-D Fe map of F07 does reveal 
some asymmetry inside 0.05r5oo, showing a high-metallicity blob 
south-west of the group core, but it is not immediately obvious 
that this is necessarily associated with an AGN outflow. To achieve 
a satisfactory fit in the outermost radial bin in our data, it was 
necessary to exclude the second brightest group galaxy, NGC 499, 
out to 6 arcmin (~ 120 kpc). Comparison to the single-T Fe profile 
of F07 shows good overall agreement, including the off-centre Fe 
peak at ~ O.OSrsoo and the decline to ^ 0.3 Zq at O.Srsoo. 

NGC 533: This group presents one of the clearest examples of an 
off-centre peak in the Fe distribution within our sample. As for 
NGC 507, which als o hosts such an Fe peak, 1.4-GHz NVSS data 
dCondon etij]|l998h reveal a central radio source. Inside 0.02r5oo 
our Fe results are consistent with the F07 value of ~ 0.7 Zq, as 
is our detection of the supersolar Fe peak at 0.04r5oo. The clear 
decline to O.I Zq at large r is not obvious from the F07 results, 
however, which show considerable scatter outside O.lrsoo. We 
speculate that this relates to our assumption of circular symmetry 
and our incomplete azimuthal coverage at large radii. The 2-D 
Fe map of F07 suggests an elongated abundance structure on 
large scales, with the ACIS-S array in our Chandra data running 
very nearly perpendicular to the major axis of this structure. This 
implies that our result includes low-metallicity gas at the largest 
radii, beyond those included in the F07 analysis, and that we are 
inevitably missing some of the higher-metallicity gas seen by F07. 

NGC 741: With a g r oup X -ray luminosity of 3 x lO"^ erg s"^ 
dOsmond & PonmanI l2004h . and a B-band magnitude differ- 



8 J. Rasmussen and T. J. Ponman 



1.0 



0.1 



1.0 r 



0.1 



1.0 
0.1 

1.00 

0.10 
0.01 

3 

^ 1.0 

] 

0.1 
1.0 



0.1 r 



1.0 



0.1 



10.0 



1.0 



0.1 



NGC 383 



Fe 



NGC 507 



NGC 533 



NGC 741 



NGC 1407 



NGC 2300 



HCG 42 



MKW 



Fe 



Fe 



Fe 



Fe 



Fe 



Fe 



Fe 



Si 



Si 



Ha 



Si 



Si 



Si 



Si 



Si 



Si 



0.001 0.010 0.100 0.001 0.010 0.100 

r/ rsoo r/ r^oo 



NGC 41 25 



Fe 



1.0 ; 



0.1 



NGC 4325 



1.0 ; 



0.1 r 



Fe 



HCG 62 

1.0 rTTTTTt— -t," 

0.1 r 



Fe 




NGC 5044 



Fe 



o 



1.0 r 



0.1 r 

; NGC 5846 



0.1 r 



Fe 



4 
+ + 



NGC 6338 



1.00 r 



0.10 r 

0.01 _ 
1.0 ; 



NGC 7619 



Fe 



Fe 



0.1 r 



0.001 0.010 0.100 0.001 0.010 0.100 



Si 



Si 



Si 



Si 



t'tl, 



Si 



Si 



Si 
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ence between the first and second brightest galaxy (the spiral 
UGC 1435) within rsoo of Ami2 = 2.6, this system meets the 
criteria for fossil groups introduced by Jones et al. ( 2003). Its 
temperature and metallicity profiles appear similar to those of 
most other systems studied here, though the Fe abundance at 



large radii is clearly lower than average. It is also a group with a 
powerful c entral radio source and large-scale radio emission (e.g., 
Ijetha et afl ^2007). although the latter is probably a tailed radio 
source associated with the nearby group member NGC 742 rather 
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than with the central elhptical itself (Jetha et al., in prep.). 

NGC 1407: This is probably a dynamically evolved group, with 
overall optical pro perties not too dissimilar to those of the fossil 
system NGC 741 ( iTrentham. Tullv & Mahdavill2006l) . Its Fe and 
Si profiles appear fairly typical for the sample. 

NGC 2300: This is the only group within our sample with a spa- 
tially resolved Fe profile outside 7325 that does not show clear evi- 
dence for a radial decline. The Si profile is also somewhat atypical 
of the sample, showing a clear rise with radius outside the core. 
The relatively low and near-constant Fe abundance, as well as the 
radial rise in S i, is in very good agreement with earlier R OSAT and 
ASCA results ( Davis et al .111 9961 : iFinoguenov et all2002h . although 
the XMM study of F06 seems to find slightly higher Fe abundances 
inside O.Srsoo than those seen here. With the central early-type 
galaxy (NGC 2300 itself) having a stella r mass ~ 10 times tha t 
of the second largest galaxy, NGC 2276 jRasmussen et al.ll2006h . 
the central galaxy could potentially have been almost solely re- 
sponsible for the chemical enrichment of the hot gas in this sys- 
tem. This is also an example of a group in which ram-pressure 
stripping of (presumab ly enriched) galactic gas is taking place 
( iRasmussen et al.ll2006l) . illustrating that an SN-driven outflow is 
not the only mechanism for enriching the hot gas in the group. 
We note that t here is also evidence of a recent merger in the cen- 
tral early-type (For bes & Thomsonll 992). with further evidenc e for 
this seen in the ICM pressure map of Finoguenov et al.l ( l2006l) . 

Given that other features of this group seem typical of our 
sample, such as the presence of a cool core and a fairly undisturbed 
X-ray morphology, the question remains how to interpret the uni- 
formly low Fe abundance and the rise in Zsi with radius. We note 
that Zpc at large radii is not atypical of the sample, nor is the ra- 
dial profile of Si/Fe. This suggests a suppression of both Fe and Si 
enrichment in the group core. In light of the evidence of a recent 
merger in the central early-type, one possible explanation is that a 
cool, low-Z system has dropped into the group core, thus diluting 
any Fe excess but preserving the cool core. 

Another possibility, potentially induced by the very low 
central temperature, is that high-Z gas has cooled out in the 
group core. Based on the measured X-ray flux, the fit t ed spe ctral 
parameters, and the gas density profile of lOavis et alj ( 1 19961) , we 
estimate a central cooling time of only '-^ 1 x 10* yr. This would 
render such cool-out a viable scenario, in particular in light of 
the absence of any strong radio source in the central early-type 
(1.4 GHz radio power P^5x 10^° W Hz"^ based on the NVSS 
flux) which could indicate ongoing AGN heating activity. In order 
to remove any Fe excess and explain the 'suppressed' Fe level in 
the core, a total Fe mass of 2.5 x 10^ Mq must have cooled out 
inside the typical radial extent of the Fe excess of r « O.OSrsoo, 
assuming an initial Fe abundance of ~ 0.7 Zq (see Section |5] 
for this choice of parameters). This is equivalent to a total gas 
mass of 2.0 x 10* M©, which does not seem an unreasonably 
large amount. In addition, cooling would also have raised the 
central entropy through the removal of low-entropy gas from the 
X-ray phase. Among the nine of our groups also included in the 
F06 and F07 samples, NGC 2300 does show the highest central 
entropy, lying 1.6a above the sample mean derived inside O.lrsoo 
(F06). One concern for the cooling scenario, howe ver, is that the 
coole d gas is likely to participate in star formation (Edwards et al.l 
l2007h . but there is no indication that the central early-type is 
unusually blue, displaying a B-K colour of 4.2. Based on the 
X-ray data alone, we conclude that the cooling scenario pro- 



vides a possible explanation for the lack of an Fe excess in this 
group, but further tests of this would require constraints on the 
amount of cold gas and young stars in the central early-type galaxy. 

HCG 42: The observed range of Zpc = 0.5 — 1.0 Zq inside 
O.lrsoo is consistent with the F07 results, as is the decline at large 
radii. However, finding a low value of ~ 0.1 Zq outside the cool 
core, we do not confirm the value of 0.2 Zq seen in the outermost 
bin of F07, certainly in part because we do not have the statistics 
required to resolve this narrow region around O.brsoo- 

MKW 4: The presenc e of an off-centre Fe peak in thi s grou p 
was already claimed bv lFukazawa. Kawano & Kawashim^ ( |2004|) . 
who furthermore noted that no radio activity or merger evidence 
has been reported for the system. Indeed, NVSS data show no 
evidence of a central radio source, nor of any larger- scale radio 
emission surrounding the central galaxy, NGC 4073 jjetha et al] 
|2007). This is in clear contrast to NGC 507 and NGC 533, where 
this Fe peak is accompanied by the presence of a central radio 
source. We confirm the high central Fe abundance of ~ 2 Zq seen 
by F07, along with the sharp decline, starting at O.OSrsoo, towards 
a value of 0.1 Zq at large radii. 

NGC 4125: This group is the lowest-temperature system in the 
sample. It shows no cool core, despite appearing reasonably relaxed 
in the X-ray (Fig.[T). Subdividing the innermost radial bin and fit- 
ting spectra for each sub-bin, each only ~ 500 pc wide, still does 
not reveal any temperature drop. As there appears to be a point-like 
X-ray source at the centre, detected above the local background 
at ~ 2a significance, one possibility is that the central tempera- 
ture rise could be due to recent AGN heating. However, although 
a radio source is present in NVSS data inside D25, this source is 
displaced from the optical and X-ray centre by 0.8 arcmin (4 kpc). 
To investigate the nature of the central X-ray source, its spectrum 
was fitted with a power-law absorbed by the Galactic value of TVh . 
This provided an unacceptable fit, even when allowing for intrinsic 
absorption beyond the Galactic value (red. > 1-7 for 9 degrees 
of freedom). In contrast, an APEC plus power-law model fits well 
(red. ~ 0.6), yielding T = 0.56 ± 0.12 and Z = O.adtoit for 
the thermal plasma, and F ~ 1.7 for the power-law, with a total, un- 
absorbed luminosity of 1 x lO''^ erg s^ ^ . The plasma parameters are 
in good agreement with those found on larger scales in the group 
core, and the power-law with a contribution from low-mass X-ray 
binaries in the central elliptical. The Chandra detection of a cen- 
tral point source might therefore simply be attributed to a strongly 
peaked surface brightn ess profile, in line with the /?OSAr finding of 
iMulchaev et alj ( |2003|) for this group. We therefore conclude that 
there is no strong X-ray evidence for ongoing AGN heating having 
raised the central ICM temperature. Further support for this inter- 
pretation is provided by the low X-ray luminosity derived for the 
very central regions, which is well below typical AGN luminosities. 

An alternative explanation for the central temperature rise 
could b e sought in the fact that the cen tral elliptical is 'dynamically 
young' jFabbiano & Schweizejl99^ . so a merger event within the 
past few Gyr could have disturbed the central regions and mixed in 
hotter gas from larger radii. The derived central cooling time is very 
short, however, ~ 1 x 10* yr as for NGC 2300, but increasing to 
more than 3 Gyr at D25 as a consequence of the sharply peaked sur - 
face brightness, and hence density, profile jMulchaev et al.ll2003h . 
As this is a factor of three above the sound crossing time within that 
region, it may suggest that the gas density distribution has largely 
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recovered from such a merger disturbance, but significant central 
cooling has yet to commence. 

Another possibility is that the situation in NGC4125 is 
somewhat similar to that seen for the fossil group NGC 6482, one 
of the coolest known fossils , with a temperature of only 0.5 keV 
outside the central 10 kpc jKhosroshahi. Jones & Ponmanll2004h . 
Despite looking relaxed in the X-ray and displaying a central 
cooling time of < 10** yr, this system also shows no cool core. A 
possible explanation is that any gas taking part in a steady-state 
cooling flow in this group is actually heated by gravitational PAV 
work as it moves inwa rds, owing to t he ver y cuspy total mass 
profile of the group. Kho sroshahi et al. I l l2004h find that, in such a 
scenario, the observed temperature profile can be explained by a 
cooling-flow mass accretion rate of just 2 Mq yr^^. The sharply 
peaked density profile seen for NGC 4125 is consistent with the 
idea that a similar mechanism could be at work in this group. 

NGC 4325: This is a group in which the extent of the metal 
excess in the core substantially exceeds the optical extent of the 
central galaxy. N o central radio source is seen in NVSS data 
( iJetha et al. 1|200'/), but there is evidence of a previous weak AGN 
outburst, with an indication that another outburst may be about 
to occur teussell. Ponman & Sanderson! 1 20071) . The observed 
near-constancy of Zyb ~ 0.7 — 0.8 Zq inside O.lrsoo agrees 
well with the F07 results, as does the very steep decline to a value 
Zpe < O.IZq atO.Srsoo- 

HCG 62: This group has b een extensively studied with ROSAT, 
ASCA Chandra , and XMM jpinoguenov & Ponmarj|l999l ; iBuotd 
l2000bl : lMoritae t al. 2006). It harbours one of the clearest example s 
of X-ray cavities in groups jVrtilek et al.l200d : lMorita et al.l2006h . 
indicating past radio activity. The older CIAO task 'MKRMF' was 
used to create spectral responses for this group, due to a lack 
of relevant calibration data for the newer task 'MKACISRMF' 
otherwise used in this study. Despite this and the complex X-ray 
morphology in th e core of this syste m, the combined Chandra 
and XMM study of iMorita et alj i200d) finds a temperature profile 
very similar to the one derived here, and the roughly solar Fe 
and Si abundances se en in the core agrees well with the results 
oflSuote (200d3) and iMorita et al.l (2006). The sharp decline in 
the Fe profile, beginning at ~ O.OSrsoo and reaching a value of 
~ 0.05 Zq at 0.6r5oo, is also in good agreement with F07. 

NGC 5044: This system has bee n studied e xtensive ly 
with both XMM and Chan dra teuote et al] l2003allbi: 
iBuote. Brighenti & Mathews! !2004h . At r « O.OSrsoo, our 
results are consistent with the solar Fe abundance found by F06. 
In our outermost radial bin, we find Zfc ~ 0.14 ± 0.04 Zq, 
confirm ing the low Fe abundance of .1-0 .15 Zp) deriv ed at large 
radii bv !Finoguenov & PoimianI ( Il999f) and !Buote et al.! ( [2004) . 

NGC 5846: The Fe map of F06 suggests a highly symmetric 
abundance structure, and supports our detection of an off-centre Fe 
peak and a decline to 0. 1 Z© at large radii. 

NGC 6338: Reaching a peak temperature of T ~ 3 keV, this is the 
hottest system in our sample. In Fig. [Tj extended emission is also 
seen from the second brightest galaxy. This has been masked out 
to a radius of 1 arcmin in all spectral fits. As for NGC 4325, the 
central metal excess in this group extends well beyond the central 
galaxy, but is here also accompanied by a central radio source. 



NGC 7619: NGC 7619 itself is one of a central pair of ellipticals of 
similar optical brightness in the group. Extended X-ray emission is 
also seen in Fig. [T] around the other pair member, NGC 7626; this 
was masked out to a radius of 2 arcmin in all spectral fits. 



4 SYSTEMATIC ERRORS 

Before proceeding to discuss the results in more detail, we address 
a number of analysis issues which could be affecting our results. 
These include the presence of Fe and Si biases in the group core, 
and the reliability of our background estimates. An immediate, al- 
though indirect, indication that such issues are not having a major 
impact on our results, is the fact, discussed below, that we detect 
the same overall radial trend in T and Z for all systems, despite 
considerable variations in both source and background flux across 
the sample. F or the temperature profi les, the agreement with the 
XMM result o f !GastaldelIo eTZH2005l) for the mean position of the 
temperature peak is encouraging. Moreover, there is good overall 
agreement with derived F e profiles for the s i x grou ps overlapping 
with the XMM sample of I pinoguenov et afl ( !2007!) . even at large 
radii where systematic uncertainties related to background subtrac- 
tion could be important. 

4.1 Geometric effects 

One issue which could call the robustness of our results into ques- 
tion is the fact that our profiles are projected ones, so all bins con- 
tain some contribution from emission further out. To test the im- 
portance of this effect, we deprojected the profiles for a few groups 
using the PROJCT model in XSPEC. For simplicity, a l-T APEC 
model was fitted to the data. The result for NGC 4325 is shown in 
Fig.jS] This group is typical of our sample in terms of the number of 
radial bins and the value of (T) ; it also shows a near-solar Si/Fe ra- 
tio throughout, so the measurements should not be biased strongly 
by fixing the relative abundances. The figure demonstrates reason- 
able overall agreement between projected and deprojected results. 
As can be seen, projection tends to slightly reduce the temperature 
gradient in the core, as would be expected. While the deprojected 
abundance profile recovers the increased Fe level in the region in- 
side ~ O.lrsoo, significant fluctuations between adjacent bins are 
evident in this region in disagreement with the projected profile. 
These are most likely artefacts related to issues of fit stability, a 
problem commonly arising for spectral deprojections involving a 
large number of bins (e.g., . Sanderson et al. 2006) . Results for less 
well-resolved systems support this suspicion and confirm the lack 
of systematic variations between projected and deprojected results. 
We therefore conclude that our results are reasonably reliable rep- 
resentations of the three-dimensional profiles. 

Another issue is the fact that we have incomplete azimuthal 
coverage at large radii. This could potentially play a role for indi- 
vidual groups in which a simple radial description of the abundance 
profile breaks down in the outskirts. As mentioned, however, we 
have attempted to select systems with a reasonably regular X-ray 
morphology, and we have also tried to accommodate any remaining 
uncertainties related to this in our error analysis. The good overall 
agreement with XMM results at large radii indicates that this effect 
is not important, at least for the groups where a comparison is pos- 
sible. In addition, any bias would be strongly suppressed when av- 
eraging the results across the group sample, which is the approach 
taken in the next Section. We therefore believe our conclusions to 
be robust against this effect. 
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Figure 5. Projected and deprojected temperature and abundance profiles for 
NGC 4325, resulting from fitting a single-temperature model to all bins. 



4.2 Observational bias in group cores 

It is well known that high-quality observations of groups in 
many cases reveal disturbed central struct ures, showing that group 
cores can be therma lly complex (e.g. iKim & Fabbianol 1 20041 : 
IO'Sullivanetalj2005h and require multi-component thermal mod- 
els to properly fit their spectra. To investigate this, we plotted the 
distribution of reduced values for all our fits, but found no clear 
trend with radius and no indication that multi-temperature models 
are generally required to describe the central regions. 

Nevertheless, particularly for systems such as NGC 383, 
where the central cool regions are not well resolved because of 
our binning criterion, the presence of multi-phase gas could lead 
to artifi cially low abundance estimates in the innermost bin tBuote 
l2000ah . and metallicities derived from l-T fits should here be re- 
garded with some caution. We tested for this effect by first radi- 
ally subdividing the central radial bin and fitting spectra extracted 
within each of these two sub-bins. Eight of the 15 groups show 
evidence for T variations between those two sub-bins. For all sys- 
tems, the abundance values within the two sub-bins were consistent 
with the full-bin value, and only one group displays evidence for Z 
variation between the two sub-bins (NGC 383, with the innermost 
sub-bin showing a lower value). However, this does not in itself 
rule out the presence of multi-phase gas within each of these two 
sub-bins. A further test therefore involved fitting two-temperature 
models to all bins within D25, as already described in Section [Z2l 
Only for five groups do we find evidence that 2-T models yield sig- 
nificantly higher abundances and provide superior fits, in some of 
the innermost radial bins (Fig. [3}. For these, the 2-T results for the 
abundances have been adopted, as plotted in Fig. [3] An example is 
NGC 507, for which only the 2-T fit recovers the high Fe abun- 



dance of 2-3 Zq reported inside D25 bv lKim & Fabbianol jiooj) 
on the basis of XMM data. 

While 2-T models certainly provide better fits in some cases, 
it is perhaps slightly surprising that they only yield significantly 
different abundances in a few of these, given the presence of cen- 
tral temperature variations in most of our groups. A possible ex- 
planation is that we are typically sampling the core regions with 
rather finely spaced annuli. The limited radial extent of each an- 
nulus may imply that mixing of temperature components may not 
be impor tant within each bin. F or example, for the very core of 
MKW 4, iFukazawa et"ai] ( |2004 also find from the same Chandra 
data that 2-r models do not produce significantly improved fits 
or higher Zfc, in ag reement with our resul ts. A different conclu- 
sion was reached bv|0' Sullivan et alj ( |2003[) based on XMM data, 
which [pukazawa et alj 1 2004 ) attribute to the much larger spectral 
extraction region used by the latter authors. Sup port for th i s inter - 
pretation is provided by the work of Humphre v & Buot3 ( l2006h . 
which shows that Chandra spectra of elliptical galaxies can gener- 
ally be described by l-T models, and that the multi-T fits required 
to describe earlier ASCA spectra of the same objects were an arte- 
fact of sampling strong temperature gradients with poorer spatial 
resolution. The fact that we generally find 2-T models to yield bet- 
ter fits and higher abundances only for groups where our innermost 
bin covers a large fraction of D25, would seem to add support for 
this interpretation. Indeed, as discussed in Paper II, for spectra cov- 
ering the full region inside O.lrsoo (a region much larger than that 
covered by the central bin in all groups), we find that two-phase 
models are in most cases statistically required. In addition, we note 
that the introduction of a F ~ 1.7 power-law to account for un- 
resolved X-ray binaries in the central early-type also mitigates the 
need for a 2-T model in some cases. 

It is finally worth noting that any remaining Fe bias would only 
strengthen the result that an Fe excess is present in the cores of most 
of our groups. Qualitatively, therefore, the only significant feature 
potentially induced by the Fe bias would be the off-centre Fe peak 
seen in some systems. However, some of our groups for which 2-T 
models provide better fits and significantly enhanced abundances 
in the central regions, such as NGC 507, still show evidence for 
such Fe and Si peaks, suggesting this feature is reasonably robust 
to such a bias. 



4.3 Background estimation 

Another important issue is the accuracy of our background subtrac- 
tion. This is mainly relevant at large radii, where the emission at all 
energies is typically dominated by the background. As mentioned, 
we use blank-sky data to extract background spectra, scaled to 
match the 10-12 keV count rates of our source data for each CCD. 
The resulting medium-to-hard background above E > 4 keV 
seems properly accounted for in all cases, as there is no statistically 
significant residual emission above this energy within the spectral 
extraction regions employed for any data set. Of potentially more 
concern is the soft X-ray background below i5 ~ 1 keV, whose 
intensity varies considerably across the sky. For 10 out of our 15 
groups, the ROSAT All-Sky Survey 7?45 (0.5-0.9 keV) background 
count rates are within 2a (< 10 per cent deviation) of the anal- 
ogous value for the relevant blank-sky background data, lending 
support to our use of the latter for background estimation. While 
most of our groups are thus not projected onto regions of signif- 
icantly enhanced or depleted soft X-ray background, there are a 
few exceptions which all show enhanced background, most notably 



Temperature and abundance profiles of groups 13 



NGC 5044 and NGC 5846 which both he close to the North Polar 
Spur. 

For most of our systems, intragroup emission covers the field 
of view (rsoo is outside the field of view in most cases), so we 
carmot in general use the outer portions of the observed fields to re- 
liably assess the level of any residual soft background. Instead, two 
independent tests were performed to assess the impact of varying 
the assumed background: We either simply rescaled the blank-sky 
background, or added a separate background model component in 
the fits to the group emission. 

In the first test, we repeated all fits for the three outermost 
radial bins of each group, with the blank-sky background normal- 
ization shifted by ±10 per cent. The fits were performed in the 
0.7-2.0 keV band only, to prevent the obvious mismatch resulting 
at higher energies from affecting the results. We found that this in 
general induced negligible changes in the best-fitting temperature, 
which remained consistent with the value derived for the 'unper- 
turbed' background in all but two cases (when increasing the back- 
ground for NGC 4125 and NGC 6338), for which the resulting fits 
were statistically unacceptable anyway. If anything, the best-fitting 
value of T is slightly lower in all cases, so there is no indication 
that the temperature decline observed in all groups at large radii 
can be ascribed to inaccurate background subtraction. We therefore 
conclude that the large-scale behaviour of the temperature profiles 
in Fig.[2]is a robust result. 

The resulting abundances are more sensitive to variations in 
the assumed background level. For three of the groups (HCG 42, 
NGC 1407, and NGC 2300) rescaling the background produced 
higher abundances in the outermost bin, inconsistent with the re- 
sult obtained using our standard background level. None of these 
cases yielded unacceptable fits, so we conclude that significantly 
higher values of Z in the outermost region probed for these systems 
are allowed, though not favoured, by the data. However, at least in 
the case of HCG 42, there is good reason to believe that the exist- 
ing abundance measurement is fairly reliable, as it compares well to 
the r esult derived independently from XMM data dFinoguenov et al.l 
l2007i) . Similar agreement is found for NGC 2300 when comp ared 
to earlier ROSAT and ASCA results (see lRasmussen et alj2006l) . As 
an additional test, we also tried varying the background normaliza- 
tion by ±20 per cent for all groups. In the four cases where the 
resulting fits did not become either significantly worse or statisti- 
cally unacceptable, the derived T and Z were consistent with the 
results obtained using the unperturbed background. 

As a second test, carried out for the five groups with signifi- 
cantly enhanced soft background, we included a spectral model to 
account for this excess background, and repeated all fits outside the 
cool core with this m odel added to that used to d escribe the group 
emission. Following IVikhlinin et al] ( |2005|) and |o: Sullivan et al.l 
( I2007 ). a fixed T = 0.18 keV, Z = Zq MEKAL background 
model was assumed, with only its normalization allowed to vary. 
Only for NGC 5846, which - along with NGC 5044 - shows the 
highest background excess, did we find significant changes in fitted 
parameters, with this (improved) fit yielding lower Fe and Si abun- 
dances, which we have consequently adopted. For NGC 5044, no 
changes were seen, possibly because of the relatively high signal- 
to-noise ratio in the outermost bins for this group. 

We confirmed that soft background excess is not an issue for 
the remainder of the sample. For example, for relatively compact 
low-surface-brightness sources such as HCG 42 and NGC 4125, 
there is no statistically significant residual emission in the 0.7- 
2 keV band in a 2 arcmin wide annulus extending to the edge of 
the field of view. For the adopted soft background model, only 



15 per cent of the total 0.5-0.9 keV flux (used to determine the 
presence of a soft excess) falls at energies above 0.7 keV. The fact 
that energies below 0.7 keV are excluded in spectral fitting should 
therefore also reduce the problem. Finally, we do not see any sys- 
tematic low-energy (~ 0.7-1 keV) residuals for any of the spectra. 
This suggests that systematic uncertainties in our background sub- 
traction is not a major issue at low energies either. 



5 COMBINING THE RESULTS 

A radial profile of some quantity may not always provide a useful 
description for an individual system due to departures from spher- 
ical symmetry. The average of many such profiles derived for dif- 
ferent groups should be much less affected though. Since our radial 
measurements for individual groups are statistically independent 
(as they are based on projected profiles), it is possible to explore 
general trends in the radial distributions of T and Z by superpos- 
ing and averaging the profiles derived for individual groups. With 
the robustness of our results well established, we can therefore pro- 
ceed to discuss these in a statistical framework o n the basis of cor- 
relation tests and orthogonal regression analyses jlsobe et aljl990l) 
performed on the combined measurements. Throughout the sec- 
tion, we have used the Kendall rank-order correlation coefficient 
r (in the range [— 1;±1], with +1 for a strong, positive correla- 
tion), and the associated significance gk of r being non-zero, to 
quantify the strength and nature of any linear relationship between 
various quantities. 

5.1 Temperature profiles 

In order to examine the overall radial behaviour of T for our sam- 
ple, we plot in Fig. |6] all temperature measurements, normalized 
to the mean temperature (T) of each group. In the plot, we have 
subdivided the data points according to whether they are inside or 
outside the T-peak located at r(rinax), resulting in 134 and 58 
measurements, respectively. For NGC 4125, which has a mono- 
tonically decreasing T(r), we have grouped the data points inside 
D25 with those at r < r(Tmax). As discussed, temperature varia- 
tions clearly exhibit a very similar overall b e haviour. Similar to the 
case for clusters (e.g. Ivikhhnin et~aLll2005l ; IPratt et al.ir2007h . and 
as discussed quantitatively below, there is a clear indication that the 
intrinsic scatter in T / (T) is larger in the central regions, inwards 
of the temperature peak. This is not surprising, as non-gravitational 
processes other than radiative cooling could be important in this 
region, thus helping to break any self-similarity otherwise induced 
by gravity. 

At first glance, Figure|6]appears to be suggestive of a bimodal- 
ity in T(r) in the core regions, but this could be largely driven by 
NGC 5044, which is better sampled than the other groups and dis- 
plays lower central values of T/ (T) than the majority. To investi- 
gate this in more detail, we tested for any evidence of bimodality 
in T/ (T) by extracting the value for each group around a specific 
overdensity radius inside the cool core, excluding the non-cool core 
system NGC 4125. For r ~ O.Olrsoo, the smallest radius shown in 
Fig. [6] a Kolmogorov-Smimov test indicates a 65 per cent probabil- 
ity that the resulting values have been drawn from a single Gaussian 
distribution with the given mean and variance. At O.OSrsoo, how- 
ever, there is some evidence for bimodality (see inset in Fig. |6ll, 
with the corresponding probability being only 12 per cent, but in- 
creasing again to 45 and 91 per cent at 0.05 and O.OTrsoo, respec- 
tively. Intriguingly, we conclude that there is some evidence for a 
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Figure 6. Temperature profiles of all groups, normalized to the mean temperature (T) of each group. Filled circles represent data points outside the position 
?'(7max) of the temperature peak (i.e. outside the cool core), empty circles those inside. Lines show results of regression fits to the data, assuming various 
functional forms of T/ (T) vs. r as described in equations j5)-0. Inset shows the histogram of T/ (T) at r O.OSrsoo for the 14 cool-core groups. 



bimodal T/ (T) distribution around r ~ O.OSrsoo, but the result is 
marginal in the sense that a unimodal Gaussian distribution cannot 
be rejected at 90 per cent confidence, in part due to the relatively 
small sample size. There is also no strong indications of bimodality 
at significantly smaller or larger radii. 

In order to investigate the potential role of AGN feedback 
in establishing this tentative bimodality, we searched for a cen- 
tral radio source in the brightest group galaxy following the pro- 
cedure outlined in ,Crosto n et al. (2005). Such a source was clearly 
identified in 11 of the groups (with the off-centre source in 
NGC4125 as an additional candidate), and with HCG42, MKW4, 
and NGC 4325 as the only clear exceptions. For these 1 1 groups, 
testing for a linear correlation between T/{T) at r « O.OSrsoo 
and the 1.4 GHz power of this central radio source yields a value 
of Kendall's r = +0.31 and an associated significance of ax = 
+ 1.3. Thus, there is marginal statistical evidence for a correlation, 
but this is largely driven by the two most powerful sources (in 
NGC 383 and NGC 741), without which any correlation is absent 
(aK = +0.4). In conclusion, there appears to be no clear connec- 
tion between this potential T{r) bimodality and the current radio 
output of a central AGN. A larger sample would be useful for a 
more decisive test of the presence of a bimodality, and for address- 
ing its possible origin. 

It is also worth noting that we are probing the cores of these 
groups down to much smaller fractions of rsoo than in typical 
cluster studies. Interestingly, excluding the non-cool core system 
NGC 4125, the ratio Tc/{T) of the temperature Tc in the innermost 
radial bin to the mean temperature (T) shows a mean and standard 
deviation of 0.58 and 0.14, respectively (0.61 and 0.17 if NGC 4125 



is included). This is higher than the typical cool-core cluster value 
of Tc/ (T) ~ 0.4 foun d, for example, for the proje cted T-profiles 
ofiVikhlinin et alj j2005i) or the de projected ones of IPiffaretti et al.l 
( l2005l) and I Sanderson et aL lllooi). It may reflect the fact that the 
temperature profiles of our groups flatten off at small radii (inside 
r « O.Olrsoo), as suggested in Section [3T1 in contrast to what 
is seen in clusters. Fig. |5] indicates that deprojecting our profiles 
would not change this conclusion, as already suggested by the in- 
ternal consistency among projected and deprojected cluster results. 

An orthogonal regression fit to the temperature profiles inside 
the cool core - excluding the innermost r < O.Olrsoo because of 
the putative flattening in this region ~ yields 

log T/{T) = (+0.21 ± 0.02) log(r/r5oo) + (0.28 ± 0.03). (3) 

The derived logarithmic slope of the profile in the core of ~ 0.2 
is in excellent agreement with the value of 0.23 ± 0.04 seen for 
groups in the radio-quiet (1.4 GHz radi o power of c entral early- 
type L < 10^^ W Hz"^) subsample of ljetha et all ( l2007h . It is, 
however, somewhat flatter than the slope of 0.34 ± 0.04 seen for 
their radio-lo ud subsample, or that of ~ 0.4 found for clusters by 
ISandersonetal.l ( l2006i) . 

Hence, not only do the temperature profiles appear to flatten 
off in the very core of these groups, their rise with radius is also 
somewhat shallower than in clusters. We investigate these results 
from a slightly different perspective in Figure|7] plotting the ratio of 
Tc to the peak temperature Tmax for each group. While this ratio is 
subject to considerable scatter at a given mean system temperature 
{T), there is a weak indication (r = —0.43, aK ~ —1.3) of a 
negative correlation between the quantities plotted (the low-lying 
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Figure 7. Ratio of core to peak temperature as a function of mean group 
temperature (T). NGC 4125, which has no cool core, is shown as an empty 
circle. Dashed line shows the typical cluster value of Tc/Tmax ~ 0.4 (e.g., 
[Sanderson et al_2006.) . 



outlier in Fig. |7] is NGC 2300; if one were to exclude this data 
point, the significance of the correlation would increase to 2.2a), 
thus confirming the notion that the cool core is less pronounced in 
groups than in clusters. The implication seems to be that, relative 
to the maximum temperature of the system, cooling in clusters can 
generally progress further than in groups, and this group-cluster 
discrepancy becomes increasingly pronounced towards lower mean 
system temperatures (T) . 

A potentially related issue is the radial extent of the cool cores. 
As mentioned in Section [3Tl the position r(rmax) of the tempera- 
ture peak varies slightly among the groups. In order to investigate 
whether this might reflect the properties of the central galaxy or the 
group itself, we plot inFig.[8]r(rn,ax) as a function of the size and 
mass (measured by D2b and Lk, respectively) of the central galaxy 
and of the mass of the group (meas ured by rspp), with JiT -band lu- 
minosities Lk taken from 2MASS ISkrutskie et alj|2006l) . Testing 
for a linear correlation between r(rmax) and the three quantities 
plotted, we find r = 0.62 and uk = 3.2, for r(rmax) vs. rsoo. The 
corresponding significances for a correlation with Lk and D25 are 
2.8 and 2.5a, respectively. Although not a strong result, it suggests 
that the central temperature structure of the group is more closely 
related to the properties of the group as a whole than to those of 
the central galaxy (any correlation with central galaxy properties 
might naturally arise because D25 and Lk are themselves corre- 
lated with rsoo at 2.5 and 2.7a significance, respectively). Such an 
interpretation would be in qualitative agreement with the conclu- 
sions of^elsdon et al. (2001) and Helsdon & Ponman (20o3) who 
found that the X-ray properties of gas surrounding central group 
galaxies are different from those of non-central galaxies and, in- 
deed, appear more closely related to the group than to the central 
galaxy itself. 

Even when excluding NGC 4125, by far the coolest group in 
our sample, a fit to the data in Fig.[8j; gives the relation 

r(rn,ax) = (0.20 ± 0.02)(r5oo/kpc) - (46.7 ± 15.1) kpc, (4) 

which, perhaps coincidentally, would predict the observed absence 
of a cool core in NGC 4125. If this relation can indeed be extended 
to rsoo ;$ 350 kpc, the implication would be that systems with 
''500 250 kpc, or, more precisely, with T < 0.3 keV from 



equation ([TJ, would generally not have cool cores. This is perhaps 
not surprising, as the central gas properties in systems of such low 
temperatures could be strongly affected by galactic feedback. Mea- 
sured relative to rsoo, Fig. [Sj; also shows that while the cool core 
is generally smaller than in clusters, it approaches the typical clus- 
ter extent at the high end of the temperature range covered by our 
sample, with group and cluster results becoming statistically indis- 
tinguishable at T > 2 keV according to equations (TJ and (O. 

Outside the cool core, the temperature profiles in Fig.[6]clearly 
decline with r, appearing to drop by almost a factor of two out to 
r5oo . A correlation test on the quantities plotted in Fig.|6]establishes 
a strong anticorrelation beyond r — r(rmax), with r = —0.70 
and aK ~ —7.7. For the data in this region, we attempted three 
different parametrizations for this decline of the T-profile. From a 
Chandra analysis of the projected temperature pro files of a sample 
of 13 massive, relaxed clusters (T-range 1-9 keV). IVikhlinin et all 
(2005) find a relation T/{T) = -0.76r/r5oo + 1.22 outside the 
cool core (assuming rsoo ~ 0.63ri8o). A similar linear fit for our 
groups yields 



T/(T) 



-0.74 ± 0.09)(r/r5oo) + (1.20 ± 0.02), 



(5) 



for a reduced Xv = 1-97. Taken at face value, this relation, shown 
as a dashed line in Fig. |6] is in excellent agreement with that of 
the Vikhlinin et al. clusters, suggesting that temperature variations 
within groups are very similar to those of clusters at large radii, 
contrary to the case in the core. The resulting fit quality is poor, 
however, so a relation of the form of equation lO is not particularly 
successful in describing our group data. The main reason is that 
such a parametrization cannot match the group data immediately 
outside the cool core without overpredicting observed temperatures 
at intermediate radii and underpredicting them further out. 

As a possible alternative, we note that cosmological hydrody- 
namical simulations invoking gravity only have suggested a 'uni- 
versal' temperature profile for massi ve (T > 3 keV) clusters, of 
the form T oc (1 + 0.75r/r5oo)"^ '^ jLoken et al.l2002l). assuming 
rsoo ~ 0.5rioo. This parametrization was shown by Loken et al] 
( |2002[) to be in good agreement with cluster results outside the core 
region, remaining a good description also when including galac- 
tic feedback and radiative cooling in the simulations. Motivated by 
this, we fitted a similar relation to our group data, yielding 

-1.45±0.14 



T/(T) = (1.29 ± 0.03)(1 + 0.75r/r5oo)' 



(6) 



for xt ~ 1-65. While resulting in a better fit than equation l|5]l, 
this relation suffers from the same shortcomings in describing our 
data, albeit to a lesser degree. Nominally, the result, represented by 
a dash-dotted line in Fig. [61 is broadly consist ent with that found 
for the simulated clusters of lLoken et al. Il l2002h . including the nor- 
malization factor of 1.33 found by the latter authors. 

Prompted by the appearance of the data in Fig. [6] a third 
parametrization of the profiles was also investigated, in which T 
is a linear function of log (r/rsoo). The fit result, 

T/(T) = (-0.51 ± 0.04) log(r/r5oo) + (0.67 ± 0.03), (7) 

is statistically better (xt ~ 1.26) than the results of equa- 
tion ^ and l[6]l. The intrinsic scatter about this relation is neg- 
ligible, whereas the corresponding scatter around equation l[3) is 
12 per cent, thus confirming the indication of Fig.[6]that intrinsic 
scatter is larger in the core regions. The derived relation of equa- 
tion l|7j is plotted as a solid line in Fig.|6] The fact that equations lO 
and ^ are good descriptions of both observed and simulated clus- 
ters, but are less representative of groups than equation Q, does 
not necessarily suggest the presence of intrinsic differences in the 
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Figure 8. Radius r(rniax) of maximum gas temperature as a function of (a) D25 of tlie central galaxy, (b) extinction-corrected K-hand luminosity Lk of the 
central galaxy, and (c) r'500 of the group. As in Fig.|7] NGC 4125 is represented by an empty circle. The solid line in Fig.[8j; (with slope = 0.2) rep resents the 
best fit to the data excluding NGC 4125, while the dashed line shows the typical cluster result of r(rniax) ~ O.lSr'soo (e.g., lvikhlimn et alj2005l) . 



large-scale temperature structure of hot gas in these systems. The 
formal agreement with cluster results in terms of best-fitting pa- 
rameters clearly indicate that the temperature profiles of cool-core 
groups at large radii are similar to those of clusters. We do not 
expect this conclusion to be affected by projection effects to an ap- 
preciable degree. Finally, we note that there is no indication that the 
possible bimodality in T / (T) at r ~ O.OSrsoo is reflected in signif- 
icantly different slopes at large r; assuming again a prescription of 
the form of equation Q outside r (Tmax), we find consistent slopes 
of -0.49 ± 0.07 and -0.54 ± 0.04 for the low- and high-r/{r) 
subsamples seen in the inset of Fig.|6l respectively. 

Summarizing, we find that cool cores in groups are less pro- 
nounced than those of clusters, in terms of both their spatial extent 
and relative temperature variation. At large radii, however, the tem- 
perature profiles of our groups decline with radius in consistency 
with results for more massive systems. 



5.2 Iron and silicon profiles 

In Fig. 121 all measurements of Zpc and Zsi have been combined, 
demonstrating more clearly the central excess in iron abundance 
and the general decrease to a value of ~ 0.1 Zq at large radii sug- 
gested by Fig. |3] For Zp^, there appears to be a clear trend with 
radius outside the cool core, and a test performed on the data points 
outside r(rmax) confirms the presence of a linear correlation be- 
tween log Zfc and log (r/rsoo) at the 5. la level (r — —0.46). 
Inside the core, the data are consistent with a constant Fe abun- 
dance (cr/f — 0.0), but this may mask a more complex behaviour 
as discussed below. 

The Si profiles are fairly similar to those of Fe in the group 
cores, showing no systematic trend with radius (ax = —0.1). Si 
also appears from Fig. |9] to show less systematic radial variation 
outside the core than Fe. Indeed, outside r(rniax) no significant 
correlation (t = —0.07, aK = —0.8) is found between log Zsi 
and log (r/rsoo), though the value of r suggests a mild overall de- 
cline with radius. The situation is qualitatively similar to that seen 
for clust ers, where Si generally shows a more uniform distribution 
than Fe jFinoguenov et al.l200(]|) . As is evident from Fig.[3] Zsi ac- 
tually rises with radius in the outskirts for some of the groups, a fact 
which clearly contributes to obscuring any statistical evidence for 
a simple linear behaviour. A similar feature has been noted in some 
relatively cool clusters but is s eemingly absent in hotter systems 
jFinoguenov et al.ll200dl2001ah . 

In addition to the considerable statistical errors on the derived 



abundance data, particularly for Si at large radii, there is clearly also 
substantial intrinsic scatter among the groups at any given radius. 
Coupled with the lack of a clear radial correlation for Si, we have 
consequently refrained from fitting regression lines to the raw data 
in Fig. |9^. In an attempt to suppress these intrinsic variations, we 
also plot in Fig. [9J5 the abundance data normalized to the mean 
value (Z). As suggested by the figure, the resulting correlation is 
slightly stronger for Zpc (t ~ —0.61, aK ~ —6.8) but is still not 
very significant for Si (aK = —1.2). The latter result is perhaps 
not surprising, given that the normalization factor, i.e. the mean 
abundance (Z), has been derived assuming all abundances fixed 
relative to solar, so the fitted value of (Z) will be driven by the 
abundance as estimated from the prominent Fe lines. 



5.3 The silicon-to-iron ratio Zsi/Zpc 

The strong radial decline of Zpc and the hint of a more gentle de- 
cline of Zsi are also reflected in an increase in Zsi/Zps at large 
radii, as illustrated in Fig.[TO] This figure further suggests the pres- 
ence of a dichotomy in the radial distribution of Zsi/Zpc, with this 
quantity being roughly constant inside r(Tmax) but rising further 
out. A correlation test confirms that within the cool core there is 
only marginal evidence for a systematic variation in this quantity 
(r = 0.08, ctk = 1.3), suggesting a mild radial decrease, with the 
data showing a mean and la dispersion pf 0.91 ±0.29 Zsi.Q/Zpc,©. 
This mean value is thus slightly lower than, but consistent with, the 
local (Solar neighbourhood) IMF and SN mixture. For the adopted 
SN yields, the value corresponds to a mixture of 70 per cent core- 
collapse SN and 30 per cent SN la by number, in excelle nt agree- 
ment with recent results for the cores of clusters (see de P laa et al.l 
I2OO7I and references therein), and consistent with the observed pro- 
portion of SN types at low redshift ( z < 0.5) determin ed from 
the Hubble and Chandra Deep Fields jPahlen et al.ll2004h . For the 
adopted SN yields, this implies that SN la are the dominant Fe con- 
tributor in group cores, being responsible for 79 ± 13 per cent of 
the Fe in this region. This agrees not only with results for cluster 
cores, as mentioned, but is also in accord with the corresponding 
fraction for the hot gas in elliptical galaxies of 66 ± 11 per cent 
( iHumphrev & Buoteil2006l) . 

Outside r(rinax), the correlation test confirms that Zsi/Zpc 
rises with radius (r = 0.21, aK ~ 2.3), but the trend is probably 
suppressed due to large statistical scatter. We note that many of the 
data points with very low Zsi/Zpc have Zsi ~ rather than high 
Zfc, but these have large fractional errors and most are consistent 
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Figure 9. (a) Zpc and Zgi, all measurements. Symbols are as in Fig. [6] (b) The same data normalized to (Z). 



with the general trend. Similar to the case for clusters and as further 
illustrated in Fig.llOb. the high central Fe abundance seen in most 
of our systems can thus be ascribed predominantly to SN la, with 
an increasing SN II contribution towards larger radii and lower Fe 
abundances. We also note that the global mean and standard de- 
viation of Zsi/Z-pc-, derived from those data points for which the 
fractional errors (at 90 per cent confidence) on both Zsi and Z^c 
are less than 100 per cent, are ~ 1.3 and 0.8 in solar units, respec- 
tively, in agreement with the (radially c onstant) ratio of 1.4 ± 0.2 
found for the T < 6 keV subsample of iTamura et alj ( |2004[) . The 
fact that our groups do not show a radially constant ratio is reflected 
in a large standard deviation around the global mean. 



5.4 Binning the profiles 

To better investigate radial trends and further suppress intrinsic 
scatter among the groups, we plot profiles of all data points in 
Fig. [TT] accumulated in radial bins of 20 measurements. In this 
representation, the temperature profiles are clearly seen in Fig. II lb 
to peak at a radius r ~ O.lrsoo, while flattening out in the group 
cores. We re-iterate that this central flattening is not observed in 
clusters, even in cases wher e the cluster cores are well resolved on 
scales below O.Olrsoo (e.g.. lSanderson et alj200 d). 

For the Fe abundance. Fig. Illb reveals roughly solar abun- 
dances in the group cores, consistent with the results found for X- 



ray bright ellipticals bv iHumphrev & Buotd {2QQ&I, five of which 
constitute the central group elliptical in our sample. The figure 
also confirms the trends discussed above, indicating that outside 
the core, Z-pc continuously declines out to the largest measured 
radii. The binned results further enable a straightforward compar- 
iso n to the corresponding Fe profile found for cool-core clusters 
bv iDe Grandi et ( |2004) . This reveals similar values of Fe en- 
richment within the central abundance excess in groups and clus- 
ters {Zyc ~ 0.8 Zq), with some indication that the excess is 
slightly broader in clusters; our binned Fe profile declines from 
its central value by a factor of two around r ~ 0.2r5o o , whereas 
the correspondin g radius for thelDe Grandi et alj j2004l) clusters is 
r ~ O.Srsoo- The lDe Grandi et al. U2004h clusters display a flatten- 
ing in the Fe profile beyond O.Srsoo- contrary to what is observed 
for our groups. Note also that the off-centre abundance peaks hinted 
at in many of the groups with well-resolved cores, and the consid- 
erable intrinsic scatter these peaks cause in this representation, are 
visible for both the binned Fe and Si profiles, even when these are 
normalized to the mean abundances. 



The binning of the abundance profiles clearly suppresses the 
considerable scatter seen in Fig. [9l resulting in a much more well- 
defined radial behaviour for both Fe and Si. We performed linear 
regression fits to these sample-averaged profiles in order to char- 
acterize the reasonably uniform behaviour outside the very group 
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Figure 10. Zgi/Zpo, all measurements, in units of the Solar ratio and as 
a function of (a) radius and (b) Zpc (with error bars on Zpo omitted for 
clarity). Symbols are as in Fig. |6] Dashed lines mark the ratios expected 
from pure SN la and SN II enrichment, as in Fig.O 

cores. Excluding the innermost data point in Fig. [TTb . orthogonal 
regression on the binned Fe profile yields 

log Zfo = (-0.66 ± 0.05) log(r-/r5oo) - (1.00 ± 0.06) (8) 

relative to the corresponding solar abundances. The Fe abundance 
at rsoo suggested by this relation of Zyc — 0.1 Zq is lower than 
the typical va lue of ~ 0.2 seen in the outskirts (at r rspp) of 
clusters (e .g.,'Finog uenov et alj200(]| : lDe Grandi & MolendilllOOll : 
IXamura et al. 2004). This is particularly true in light of the fact 
that w e have adopted the abundance table of iGrevesse & Sauvall 
(Il998l). for which the Sola r Fe abundance is only 68 per cent of the 
lAnders & Grevess3 ( Il989h value typically used in previous cluster 
studies. The Fe abundance in those previous studies must there- 
fore be multiplied by 1.48 for a direct comparison to our results 
(the Si abundance, on the other hand, remains unaffected to within 
1 per cent by this change in abundance table). 

Figlllt suggests that Zsi also declines outside the cool core, 
despite the fact that there is only marginal statistical evidence for 
a radial decline of Zsi outside r(Tnmx) in the unbinned data. As 
already indicated by the dichotomy in the Si/Fe ratio (Fig.llOt. the 
reason for this apparent inconsistency is that while Si does decline 
immediately outside the cool core, it flattens out or even rises at 
the largest radii probed, thus suppressing any statistical evidence 



for a uniformly declining profile outside the core. Beyond 0.4r5oo, 
roughly corresponding to the range covered by the outermost data 
point in Fig. Illb . the unbinned Si data show no clear correlation 
with r {ok ~ —0.7). To capture this behaviour, we performed 
a regression fit to all but the innermost data point in Figll lb. and 
assumed Zsi (r) to remain constant beyond the outermost data point 
at r « 0.5r5oo, yielding 

i„„ 7, ^ / "ilog(r/r5oo) +c, r < O.Srsoo ™ 
® ^' I 0.25 ± 0.05 Zq , r> 0.5r5oo ' 

with m = -0.44 ± 0.07 and c = -0.74 ± 0.06. The value of 
0.25 Zq at r > O.Srsoo was chosen to ensure continuity, with its 
errors based on the fractional errors on the outermost data point. 
An Si abundance of ~ 0.2-0.25 Zq at rs pp is also in good agree - 
ment with earlier A5CA results for groups jFinoguenov et al.l2002h . 
whereas cluster s tend to s how slightly higher values of ~ 0.3- 
0.5 Zq (Fino guenov et al.l i2000. 2001a) at large radii. This sug- 
gests a higher content of both Fe and Si in the outskirts of more 
massive systems. However, a direct comparison is somewhat im- 
peded by the large statistical uncertainties on Zsi in groups, as well 
as the fact that Zsi{r) does not show a uniform behaviour at large 
r across our sample. 

For the Zsi/Zpa ratio, a subsolar value is seen in the inner- 
most bin in Fig. Ill F. with the ratio rising to ~ 1 within O.lrsoo 
and further to > 2 well outside the central galaxy. As discussed 
in Section|4] we believe that the result for the outermost bin is not 
seriously affected by uncertainties in our background subtraction. 
The result demonstrates more clearly what was already indicated 
by Fig. [To] that SN la provide a dominant, but not exclusive, con- 
tribution to the enrichment within the cool core and hence inside 
the optical extent of the central galaxy. Well outside the group core 
the dominant contribution comes from SN II, and at large radii the 
Zsi/Zpe ratio rises to become consistent with the value of ~ 2.6 
expected for pure SN II enrichment. Corrected for differences in 
adopted abundance tables, th e value of the outermost bin in Fig. II If 
is at least as high those of the lFinoguenov et al] l l200(]h clusters. Our 
results, therefore, do not support the claim of the latter authors that 
groups on average have lower Si/Fe ratios at large radii than clus- 
ters. The results further indicate that, on average, both SN la and 
SN II enrichment is required at all radii within at least r ~ 0.5r5oo. 
While this relies on the assumption that derived abundances in the 
central region are not seriously affected by the Fe and Si biases, 
similar central Zsi jZ-pc ratios of > 1 have also been found in ASCA 
and XMM observat ions of hotter clusters ( Finoguenov et alJl200Cl : 
iTamura et alj2004h . although it should be mentioned that the cores 
of these clusters are generally less well resolved than those of our 
groups. 

Summarizing, we find that the distributions of both Fe and Si 
in the unbinned data (Fig.|9]l appear roughly uniform in the group 
cores, but the averaged profiles clearly confirm the evidence for off- 
centre peaks in some of the groups. The enrichment in the cores 
contain a significant contribution from SN la, probably hosted by 
the central, bright galaxy, an issue to which we will return in Pa- 
per II. On average, SN II dominate by number at all radii, in agree- 
ment with the results obtained by iFinoguenov & PonmanI ( 1 19991) 
for a sample of three groups, and can on their own explain the ob- 
served abundances well outside the group core. We believe this lat- 
ter conclusion to be reasonably robust with respect to any remain- 
ing background issues. A high Si/Fe ratio consistent with, or even 
exceeding, the expectation from our adopted SN II yields is seen 
in the outermost bin of all of our groups, including groups such 
as NGC 383, HCG 62, and NGC 5044, for which the data feature 
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Figure 11. Profiles for all groups, binned into radial bins of Af = 20 data points. EiTors represent the standard deviation in each bin divided by vJV. Dashed 
lines in (a)-(c) mark the best-fitting regression lines across the radial ranges shown, i.e. equations (5), (Sj, and (9j. 



high signal-to-noise ratios e ven at large radii. Suppor t for this re- 
sult comes from the work of iFinoguenov et al] ( l2000f ). in which a 
rise in Zsi/^Fc to supersolar values at large radii is also seen for 
the groups in their sample. The trend of Si/Fe increasing with ra- 
dius has also been reproduced in cosmological simulations in which 
the galactic feedback implem entation favours the ejection of SN II 
products over those of SN la dRomeo et alj|2006l) . 

We note that the interpretation that SN II alone can explain the 
observed abundance pattern at large radii is reasonably robust to the 
choice of theoretical SN II yields available in the literature. For ex- 
ample, for a st andard Salpeter IMF , none of the SN II model yields 
considered bv I Gibson et alj ^19971) would predict Zsi/^pc > 2.6. 
Hence our conclusion regarding the SN II contribution is only com- 
promised if invoking a more 'top-heavy' IMF, with x < 1.0 in 
equation In connection with this, it is worth pointing out, on the 
basis of the outermost three data points of Fig. [TIT , that the slope 
of Z{r) for at least one of the two elements must change beyond 
'"500 if the observed Si/Fe ratio is to remain consistent with the IMF 
and SN model yields adopted here. Deep, high-S/N X-ray observa- 
tions of group outskirts could shed light on this issue and thus help 
constrain theoretical SN yields, but such a study could well be be- 
yond the capabilities of current X-ray telescopes for all reasonable 
exposure times. With their higher X-ray surface brightness at large 
radii, cool to intermediate-temperature clusters (T ~ 2 — 4 keV) 
remain more attractive targets for this purpose, also because they 
show promi nent emission lines from a wid er range of relevant ele- 
ments (e.g., IFinoguenov et al.ll2000ll2001ij) . 



6 SUMMARY AND CONCLUSIONS 

We have conducted a homogeneous analysis of a sample of 15 X- 
ray bright galaxy groups observed by Chandra, selected for their 
relatively undisturbed X-ray morphology and good photon statis- 
tics. Correlation tests and regression fits have been performed on 
the combined radial temperature and abundance profiles and de- 
rived quantities, both to investigate the nature of any relationship 



between various quantities, and to enable straightforward compari- 
son to results of numerical simulations of group formation and evo- 
lution. 

Our derived radial temperature profiles show remarkable sim- 
ilarity when scaled to the mean system temperature (T) derived 
outside the central region, and reveal that all but one group show 
clear evidence for a cool core. We find that these cores are smaller 
than those seen in clusters, both in absolute physical extent and rel- 
ative to a fixed overdensity radius, and that their sizes appear more 
closely related to the properties of the group as a whole than to 
those of the central early-type galaxy present in all our groups. The 
absolute and relative cool core sizes approach those of clusters at 
the high-temperature end of our sample. The temperature profiles 
in the core are also shallower than those seen in clusters, equiva- 
lent to a less pronounced temperature decline in the core relative 
to the peak temperature displayed by the system. This may indi- 
cate that radiative cooling can generally progress further in more 
massive systems. The flatter group profiles inwards of the tempera- 
ture peak relative to those of more massive systems certainly point 
to clear differences in the central ICM heating history among sys- 
tems of varying mass. We find tentative evidence for a bimodality in 
T/(T) inside the cool core, but this is confined to a narrow interval 
in radius and is not mirrored by differing slopes in the temperature 
profiles at large radii, which we find to be similar to those found for 
both observed and simulated clusters. 

The Fe and Si abundance profiles show a number of distinct 
features common to most of the groups. Similar to the case for cool- 
core clusters, a central enhancement in Fe abundance is clearly seen 
in all but two groups, one of which is NGC 4125, the only non- 
cool core group in our sample. This Fe excess can extend well be- 
yond the optical extent of the central early-type galaxy present in 
all groups, and is in many cases accompanied by a similar feature in 
the Si distribution. In many of the groups with a well-resolved core, 
evidence for an off-centre abundance peak is also observed inside 
the central galaxy, with a decline in abundance at the very centre, 
similar to results for some well-resolved cluster cores. Outside the 
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central excess, the Fe profiles generally show a steep decline, ap- 
proaching a value of Zpe ~ 0.1 solar at rsoo- This is lower than 
the corresponding typical cluster value by a factor of two, suggest- 
ing fundamental differences in enrichment history for the outskirts 
of these low-mass systems compared to more massive clusters. Si 
shows similar overall features but declines less steeply with radius, 
and flattens out or even rises again at large radii for the majority of 
the groups. On average. Si reaches a value of Zsi ~ 0.25 — 0.3 so- 
lar at rsoo, again slightly below typical cluster values. 

For the adopted SN yields, the ratio of Si to Fe abundance in 
the groups reveals a clear predominance of SN la enrichment in 
the central regions, with 79 ± 13 per cent of the Fe originating in 
SN la. This is in line with results for both hot gas in elliptical galax- 
ies and the cores of massive clusters, suggesting a very similar en- 
richment istory for the central regions in such systems across three 
orders of magnitude in total system mass. The SN II contribution in 
our groups increases with radius, however, and at the largest radii 
probed, or, equivalently, at the lowest Fe abundances, the Si/Fe ra- 
tio is close to, and nearly always consistent with, the expectation 
for pure SN II enrichment in all groups. We find no indication that 
Si/Fe ratios at large radii are significantly different from the corre- 
sponding values in clusters. We note, however, t hat the existence 
of two separate SN la populations, as suggested bv lMannucci et al.] 
( l2006h . could complicate the interpretation of these results. 

The abundance distribution in the intragroup medium is ex- 
pected to be tightly linked to the star formation history in groups. In 
a forthcoming companion paper, we decompose the derived abun- 
dance profiles into contributions from the two major supernova 
types to explore this issue, along with the implications of our results 
for supernova and ICM energetics, the role of AGN and mergers in 
redistributing and reheating the enriched gas, and more generally 
for the enrichment history of the groups. 
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